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Abstract 



In the type-I seesaw mechanism, both the light Majorana neutrinos (z/ l5 z^ 2 , ^3) 
and the heavy Majorana neutrinos (N 1 , • • • , N n ) can mediate the neutrinoless 
double-beta (Qv/3/3) decay. We point out that the contribution of v i to this 
0vf3f3 process is also dependent on the masses M k and the mixing parameters 
R ek of N k as a direct consequence of the exact seesaw relation, and the effective 
mass term of v i is in most cases dominant over that of We obtain a new 
bound I S-Refc-MjJ < 0.23 eV (or < 0.85 eV as a more conservative limit) at 
the 2a level, which is much stronger than \ R 2 e kM k l \ < 5 x 10~ 8 GeV -1 
used in some literature, from current experimental constraints on the 0u(3(3 
decay. Taking the minimal type-I seesaw scenario for example, we illustrate 
the possibility of determining or constraining two heavy Majorana neutrino 
masses by using more accurate low-energy data on lepton number violation 
and non-unitarity of neutrino mixing. 
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I. INTRODUCTION 



Almost undebatable evidence for finite neutrino masses and large neutrino mixing angles 
has recently been achieved from solar, atmospheric, reactor and accelerator neutrino oscilla- 
tion experiments [1-4] . This exciting breakthrough opens a new window to physics beyond 
the standard model (SM), since the SM itself only contains three massless neutrinos (i.e., 
z/ e , v and u T , corresponding to the mass eigenstates v x , v 2 and u 3 ). The simplest way to 
generate non-zero but tiny neutrino masses m i for u i is to extend the SM by introducing at 
least two right-handed neutrinos and allowing lepton number violation. In this well-known 
(type-I) seesaw mechanism [5], the SU(2) L x £7(1) y gauge-invariant mass terms of charged 
leptons and neutrinos are given by 

"Anass = kY t HE R + T h Y v HN K + l -NlM K N K + h.c. , (1) 

where H = i<j 2 H*, /l denotes the left-handed lepton doublet, and M R is the mass matrix 
of right-handed neutrinos. After the spontaneous gauge symmetry breaking, we arrive at 
the charged-lepton mass matrix M l = Y t v and the Dirac neutrino mass matrix M D = Y u v, 
where v ~ 174 GeV is the vacuum expectation value of the neutral component of the Higgs 
doublet H and characterizes the Fermi scale of weak interactions. The mass scale of M R 
(or equivalently the seesaw scale A ss ) is crucial, because it is relevant to whether the seesaw 
mechanism itself is theoretically natural and experimentally testable. Although A ss <C v 
is not impossible [6], it is in general expected that A ss should be much higher than the 
Fermi scale. In particular, the conventional seesaw mechanism works at a scale which is not 
far away from the scale of grand unified theories. Driven by the upcoming running of the 
Large hadron Collider (LHC), more and more attention has been paid to the TeV scale at 
which the unnatural gauge hierarchy problem of the SM may be solved or softened by new 
physics. If the TeV scale is really a fundamental scale, then we are reasonably motivated to 
speculate that possible new physics existing at this scale and responsible for the electroweak 
symmetry breaking might also be responsible for the origin of neutrino masses. In this 
sense, it is meaningful to investigate the TeV seesaw mechanism and balance its theoretical 
naturalness and experimental testability at the energy frontier set by the LHC [7]. 

A direct test of the type-I seesaw mechanism demands the experimental discovery of 
heavy Majorana neutrinos N k (for k = l,---,n) at the LHC, but two prerequisites must 
be satisfied: (a) their masses M k must be of (9(1) TeV or smaller; and (b) their couplings 
to charged leptons R ak (for a = e,fj,,r and k = l,---,n) must not be too small. The 
strongest bound on M k and R ek comes from the non-observation of the neutrinoless double- 
beta (Ovfifi) decay [8], as N k can mediate this lepton-number- violating process. Current 
experimental lower limit on the half-lifetime of the 0v(3{3 decay is usually translated into 



n p2 



< 5 x 1(T 8 GeV" 1 (2) 



in some literature [9]. In obtaining Eq. (2), one has ignored the contribution of three light 
Majorana neutrinos v i (for i = 1,2, 3) to the 0u(3(3 decay. 
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The first purpose of this paper is to point out that the constraint in Eq. (2) is not 
always useful for the type-I seesaw mechanism either at a super high-energy scale or at the 
electroweak or TeV scale. The reason is simply that the contribution of u i to the Ovfifi 
decay is in most cases dominant over the contribution of N k to the same process, leading to 
a much stronger bound on M k and R ek through the exact seesaw relation: 

< 0.23 eV (or < 0.85 eV) (3) 

at the 2cr level, which is equivalent to (m) ee < 0.23 eV (or < 0.85 eV as a more conservative 
bound) [10,11] for the effective mass of the Ov/3/3 decay mediated by v i . The second purpose 
of this paper is to look at whether the future measurements of lepton number violation and 
non-unitarity of neutrino flavor mixing are possible to shed light on M k . Taking the minimal 
type-I seesaw scenario [12] for example, we shall illustrate the possibility of determining or 
constraining two heavy Majorana neutrino masses by using more accurate low-energy data 
on the OvfiP decay and non-unitary neutrino mixing and CP violation. 



£ 

fc=i 



R 2 ek M k 



II. STRONGER BOUND ON THE OuP/3 DECAY 



After the spontaneous gauge symmetry breaking (i.e., S77(2) L x U(1) Y 
mass terms in Eq. (1) turn out to be 



-C! 



E L M l E R + -(v L N£) 



M, 

ME 



D 



(K 
U R 



+ h.c. , 



U(1) C J, the 



(4) 



where E and u L stand respectively for the column vectors of (e, //, r) and (is e , v , ^ r ) L . With- 
out loss of generality, one may take M l = Diag{m e , m , m T }. The overall (3 + n) x (3 + n) 
neutrino mass matrix in Eq. (4) can be diagonalized by a unitary transformation; i.e., 



M, 



ME 



D 



R 



V R 

s u 



M u 






M, 



N 



V R 

s u 



(5) 



where M u = Diag{m 1? m 2 , m 3 } and M N = Diag{M 1? • • • , M 3 }. After this diagonalization, 
the flavor states of three light neutrinos 
mass states of light and heavy neutrinos 



can be expressed in terms of the (3 + n) 
3 and Ni, ■ ■ ■ , N n ), and thus the standard 



charged- current interactions between v a and a (for a = e, //, r) can be written as 
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W~ + h.c. 

fjj 



(6) 



in the basis of mass states. So V is just the neutrino mixing matrix responsible for neutrino 
oscillations, while R describes the strength of charged-current interactions between (e, /i, r) 
and (N 1 , • • • , N n ). V and R are correlated with each other through VV^-\-RRJ = 1. Hence V 
itself is not exactly unitary in the type-I seesaw mechanism and its deviation from unitarity 
is simply characterized by non- vanishing R [13]. 
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Note that V and R are also correlated with each other through the exact seesaw relation 



VM V V 2 



RM N R 



T 



o, 



(7) 



which can directly be derived from Eq. (5). Taking the (ee) -elements for both terms on the 
left-hand side of Eq. (7), we immediately arrive at 



ei m i 



(rm n r t ) 



E* 

fc=i 



2 ek M k 



(8) 



This simple but interesting result implies that the effective mass of three light Majorana 
neutrinos in the 0v(3{3 decay is directly associated with the masses, mixing angles and CP- 
violating phases of heavy Majorana neutrinos in the type-I seesaw mechanism: 



( m )ee = 



rrii 



E RlkM k 



k=l 



(9) 



Note that both light Majorana neutrinos v i and heavy Majorana neutrinos N k can mediate 
the Ov/3/3 decay, as shown in Fig. 1. When the contribution of N k is least suppressed [14], 
the overall decay width of the Qv(3(3 process in the type-I seesaw scenario can approximately 
be expressed as 1 



i=l 
n 

E Rl k M k 



n r->2 

E TT M aF( A , M k 



k=i 



k=i M k 
Ml 



1 + 



Ml 



(10) 



where A is the atomic number, J- {A, M k ) ~ 0.1 depending mildly on the decaying nucleus, 
and M A ~ 900 MeV [14]. Given M k > 10 2 GeV, the second term in the square brackets of 
Eq. (10) turns out to be < 8.1 x 10~ 6 . Hence this term is negligible in most cases, unless the 
contribution of v i is vanishing or vanishingly small due to a contrived cancellation among 
three different V^ i m i terms (or equivalently, among n different R\ k M k terms), which is in 
principle not impossible. Let us consider two limits in which the contributions of light and 
heavy Majorana neutrinos to T Qu/3l3 are decoupled. 



In the limit of ^Ik^k 1 F(A, M k ) — > 0, which is almost a realistic case, Eq. (10) is 



k=i 



directly simplified to 



Current calculations of the nuclear matrix elements of the 0^/3/3 decay are quite uncertain due 
to our poor knowledge of the nucleon wave functions [8]. For point-like nucleons, a Yukawa-type 
potential has been used to estimate the nuclear matrix elements [14]. The relevant result remains 
true even if the finite size of the nucleons is taken into account [14]. We thank the referee for 
calling our attention to these points. 
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n 



Qu/3/3 



oc (m) 



ee 



E ^e A 



(11) 



fc=l 



Current experimental data on the 0is[3[3 decay yield an upper bound on this effective 
mass term, (m) ee < 0.23 eV at the 2a level [11], which has extensively been used 
to constrain the masses, flavor mixing angles and Majorana CP-violating phases of 
three light neutrinos in the unitary limit of V. However, it should be kept in mind 
that this upper bound corresponds to some "favorable" values of the relevant nuclear 
matrix elements [8]. If their "unfavorable" values are used, one may also arrive at 
(m) ee < 0.85 eV at the 2a level [11]. 



as a rough approximation. Imposing the bound (m)' ee < 0.23 eV and inputting M 
900 MeV and T(A, M k ) ~ 0.1 [14], for example, we obtain 



which is a bit stronger than the upper bound shown in Eq. (2). If the more conservative 
bound (m)' ee < 0.85 eV is taken, one will arrive at E-R^MjT 1 < 1-° x 1°~ 8 GeV" 1 , 
much closer to the result given in Eq. (2). Such rough bounds have been used by a 
number of authors in their preliminary studies of possible collider signatures of heavy 
Majorana neutrinos [9]. 

Note again that we have ignored the mild dependence of T{A, M k ) on the decaying nuclei 
in the above discussions. Otherwise, different 0u(3(3 decays should be separately analyzed. 

Below Eq. (10), we have pointed out that the contribution of three light Majorana 
neutrinos v { to T 0u/3/3 is dominant in most cases. This observation is especially true for the 
conventional type-I seesaw mechanism with superhigh M k (e.g., max (M k ) ~ 10 15 GeV) and 
extremely small R ak (e.g., \R ak \ ~ 10~ 13 ) [15]. When the seesaw mechanism is realized 
at the electroweak or TeV scale to generate experimentally accessible signatures of heavy 
Majorana neutrinos N k at the LHC, however, one usually has to require max (M k ) < 0(1) 
TeV and \R ak \ ^ 10~ 3 up to 0(0.1) [7], which imply a terrible cancellation in every term of 
Eq. (8) so as to give rise to tiny masses of v { . Because such a terrible cancellation in (m) ee 
does not necessarily mean the same cancellation in (m)' ee , it is possible to get (m)' ee ^> (m) ee 
special already discussed in Eqs. (12) and (13). But the situation might become 

quite subtle if the masses of heavy Majorana neutrinos are degenerate [16]. Since the function 
J-"(A, M k ) depends both on the atomic number A and the heavy Majorana neutrino masses 
M k , {m)' ee might be exceedingly small for one decaying nucleus in the (m) ee limit but 



n 



• In the limit of ^ Rlk^k ~^ 0> which is rather contrived, Eq. (10) can be simplified to 

k=i 





-i 



(13) 
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not for another in the same limit. A careful analysis of the relative magnitudes of (m) ee and 
{m)' ee for different Qv(3(3 decays is nevertheless beyond the scope of the present paper and 
will be done elsewhere. 

For n = 3, R can be parametrized in terms of nine rotation angles 9^ and nine phase 
angles 5^ (for i — 1, 2, 3 and j = 4, 5, 6) [13]. In this representation, 

(m) ee = |M 1 ^ 4 44e- 2 ^+M 2 44e- 2 ^+M 3S ? 6 e- 2 ^e 

« |M 1 s 2 4 + M 2 s 2 5 e 2i(<5 i4- <5 i5) + M 3 s 2 6 e 2i(<5 i4-' 5 i6)| ? (14) 

where = sin 9^ and = cos 9^. The approximation made in Eq. (14) is very reasonable 
because \RB)\ is at most of (9(1CT 2 ) [17] and thus all the mixing angles of R must be very 
small (at most at the 0(0.1) level). Given (m) ee < 1 eV and M 1 « M 2 « M 3 ~ u or 0(1) 
TeV, for instance, the constraint in Eq. (14) implies that two phase differences S 14 — 5 15 and 
S 1A — 8 l6 should be very close to ±7r/2 in order to assure significant cancellations among three 
terms. Fixing M i ~ 10 2 GeV (for % — 1,2,3) as an example, we find that the level of fine- 
tuning is at least of C(10~ 9 ) with s\ ~ 10" 2 or C(10~ 7 ) with s\ ~ 1CT 4 (for j = 4,5,6). 
Such unnatural cancellations seem to be unavoidable in the type-I seesaw models at the 
electroweak or TeV scale, unless the relevant mixing angles are extremely small. Current 
experimental data can only provide us with a rough bound s\ A + s 2 5 + s\ & ^ 1 . 1 x 1CT 2 [17], 
unfortunately. In a specific type-I seesaw model with the fine-tuning conditions (M D ) 3i oc 
(M D ) 2l oc (M D ) U (for i = 1,2,3) and (M D ) 2 1 /M 1 + (M D ) 2 2 /M 2 + (M D ) 2 3 /M 3 = [18], one 
may obtain M v = —M D M R M^ = 0, which in turn leads to {m) ee = 0. Then non-zero but 
tiny M v and (m) ee can be achieved by introducing a small perturbation to the texture of 
M D . If the magnitudes of M k are too big or those of R ak are too tiny, of course, there will 
be no hope to produce and detect heavy Majorana neutrinos and test the seeasw mechanism 
at the LHC [19]. 



III. THE MINIMAL SEESAW SCENARIO 

The exact seesaw relation in Eq. (7) allows us to determine M k in terms of m i and the 
mixing parameters of V and R. To illustrate this point, let us focus on the minimal type-I 
seesaw scenario which contains only two heavy Majorana neutrinos [12]. In this simpler 
case, it is easy to obtain two real and linear equations of M 1 and M 2 from Eq. (8) 2 : 

2 3 

£Rei? 2 fc M fc = - Re Ki m i » 
k=i i=i 

j2lmR 2 ek M k = -^ImVjm,. (15) 

k=l i=l 



2 Similar equations of M 1 and M 2 can also be obtained from the exact seesaw relation in Eq. (7). 
A detailed analytical calculation and numerical analysis of the dependence and consequences of 
such equations will be presented elsewhere [20]. 
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Note that either m l = or m 3 = must hold in the minimal type-I seesaw model [12], and 
thus the non-vanishing neutrino masses can be determined from current experimental data 
on two independent neutrino mass-squared differences Am^ = m\—m\ and Am^ = ml— ml 
corresponding to solar and atmospheric neutrino oscillations. After a simple calculation, we 
arrive at 



M ± = -f 



3 3 

Rei? 2 2 £ ImV^m, - lmR 2 e2 £ ReV'm, 

i=l i=l 



ReR z el lmRi 2 - lmR z el ReR z e2 

ReR 2 el £ Im^X - Imi& £ R^m, 

M 2 = ^ ^ „ 1=1 „ . (16) 

2 Rei^Imi^ - Imi^Rei^ V ; 

Using the exact and convenient parametrization of V and R advocated in Ref. [13], we have 

Kl = C 12 C 13 C 14 C 15 ' K2 = S 12 C 13 C 14 C 15 e 12 ' ^e3 = S 13 C 14 C 15 e 13 ' 

R e i = s u c 15 e~^ , R e2 = s 15 e'^ (17) 

for the minimal type-I seesaw scenario, where c u = cos 6*^ and s u = sin 6*^ (for i — 2, • • • , 5). 
There are at least two phenomenological merits of this parametrization for our present 
discussions: (1) it can automatically reproduce the standard (unitary) parametrization of 
the light neutrino mixing matrix [21] when the non- unitary mixing angles of R are switched 
off; and (2) it can lead to a very simple result of (m) ee , which is equal to the standard 
(unitary) expression of (m) ee multiplied by a factor c 2 4 c 2 5 [13]. Substituting Eq. (17) into 
Eq. (16), we obtain the explicit results of M 1 and M 2 for two different patterns of the light 
neutrino mass spectrum. 



A. Normal hierarchy (m 1 = 0) 



In this case, it is straightforward to obtain m 2 = y Am 21 m 8.8 x 10 3 eV and m 3 = 

Am 2l + |Am 2 2 | 5.0 x 10~ 2 eV, where we have typically input Am 21 = 7.7 x 10~ 5 eV 2 
and |Am 2 2 | = 2.4 x 10~ 3 eV 2 [11]. The expressions of M 1 , M 2 and (m) ee are given by 

m 2 s\ 2 c 2 l3 sin (0 2 + 0) + m 3 s\ 3 sin (0 2 - 0) _ c 2 ^ 

sin (0 2 - X ) s 2 4 ' 

m 2 s\ 2 c 2 13 sin (0 X + 0) + m 3 s 2 13 sin - 0) c 2 4 c 2 5 
sin (0 2 - X ) s? 5 

c i4 c i5 \/ m 2Si2Ci3 + m isi3 + 2m 2 m 3 s 2 12 cl 3 s 2 13 cos 20 , (19) 

where = <5 13 — 5 12 , 1 = 25 14 — (5 12 + <5 13 ) and 2 = 25 15 — (5 12 + 5 13 ). Note that M x > 
and M 2 > require that three CP-violating phases 0, 1 and 2 should not all be vanishing; 
instead, they must take proper and nontrivial values. 



M 1 



Mo 



together with 



(m) e 
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The above results can be simplified in the limit of s 13 — > 0, which actually has no conflict 
with current experimental data: 



(m) ee sin (0 2 + 0) 



Mj = 



M 2 = +^ ■ g^fi-^ > (20) 



sf 4 cj 5 sin ^ 2 

Wee _ Sin (0! 

s?5 sin (02-01^ 



and (m) ee = s\ 2 c^ A c\ b m 2 . We have the following observations: (a) the magnitude of {m) ee 
is of (9(10 -3 ) eV [22], and thus it is experimentally inaccessible in the near future; (b) the 
signs of sin(0 2 + 0) and sin(0 2 — X ) must be different, while the signs of sin(0 1 + 0) and 
sin(0 2 — X ) must be the same; and (c) it will be possible, at least in principle, to obtain 
the lower bounds on M l and M 2 if some constraints on the non-unitary mixing angles and 
CP-violating phases become available at low energies. To achieve M 1 ~ M 2 ^ 10 2 GeV, 
for example, Eq. (20) implies that s 2 4 c 2 5 sin(0 2 — X ) ~ s 2 5 sin(0 2 — : ) P$ 10~ 14 must hold. 
Current experimental bounds can only provide us with s 2 4 + s 2 5 < 1.1 x 10~ 2 [17]. Hence 
there is a long way to go before we can constrain s 14 , s 15 and even the relevant CP-violating 
phases to a much better degree of accuracy. 



B. Inverted hierarchy (m 3 = 0) 



In this case, it is easy to obtain m 1 = y |Am 32 | — Am 21 ~ 4.8 x 10 2 eV and m 2 = 

^\Am 2 32 \ « 4.9 x 10~ 2 eV, where Am 2 21 = 7.7 x lO" 5 eV 2 and |Am 2 2 | = 2.4 x 10~ 3 eV 2 [11] 
have typically been input. After a straightforward calculation, the explicit expressions of 
M 1) M 2 and (m) ee are given as follows: 

mjc|g sin (0 2 - 0Q + m 2 s\ 2 sin (0 2 + 0Q _ c 2 3 c 2 4 

sin (02 - 0i) s 2 4 
mjclg sin (0^ - ; ) + m 2 s 2 2 sin (0' x + ; ) c 2 3 c 2 4 c 2 5 

sm(0 2 -0i) S ? 5 ' 1 J 



= c i3 c i4 c i5V m i c i2 + m i si 2 + 1m x m 2 c\ 2 s\ 2 COS 20' , (22) 

where 0' = —5 12 , 0' x = 25 14 — S 12 and 0' 2 = 25 15 — 5 12 . The requirement of M : > and M 2 > 
implies that three CP-violating phases 0', 0' : and 0' 2 cannot all be vanishing; instead, they 
must take proper and nontrivial values. 

One can see that M 1 , M 2 and (m) ee are independent of the small mixing angle 9 13 . 
Given 9 12 ~ 34° [11], for example, the maximal value of (m) ee in Eq. (22) is achievable at 
0' ~ 0; i.e., {m) ee can maximally amount to (m^^ + tti 2 s1 2 ) ~ 4.9 x 10~ 2 eV [22], which is 
experimentally accessible in the near future. Again, we may optimistically argue that it is 
in principle possible to obtain the lower bounds on M x and M 2 if some constraints on the 
non-unitary mixing angles and CP-violating phases become available at low energies. 



M x = 
M 2 = 

( m )e, 
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It is worth remarking that the non-unitarity of V is signified by non- vanishing R, whose 
mixing angles and CP-violating phases are quite possible to have some nontrivial observable 
effects. For example, an appreciable CP-violating asymmetry up to the percent level is 
expected to show up between v — > u T and V V T oscillations in some medium- or long- 
baseline experiments at a neutrino factory [23], just as a consequence of non- vanishing R. 
A neutrino telescope could also be a useful tool to probe the non-unitary effect in ultrahigh- 
energy cosmic neutrino oscillations [24]. 

IV. SUMMARY 

We have carefully examined the contributions of both light Majorana neutrinos v i with 
masses m i (for % = 1,2,3) and heavy Majorana neutrinos N k with masses M k (for k = 
l,---,n) to the OvfiP decay in the type-I seesaw mechanism, in which the light neutrino 
mixing matrix V is non-unitary due to the non- vanishing coupling matrix R between N k 
and charged leptons. The exact seesaw relation allows us to establish a straightforward 
relationship between (m i: V ai ) and (M k) R ak ). We have pointed out that the constraint 
\Y,R 2 e kM k ~ l \ < 5 x 10~ 8 GeV -1 used in some literature is in most cases too loose for a 
type-I seesaw mechanism either at a superhigh-energy scale or at the electroweak or TeV 
scale, because the contribution of v i to the Ou/3/3 decay is in most cases dominant over the 
contribution of N k to the same process. Such an observation leads us to a much stronger 
bound on M k and R ek ; i.e., | J2R 2 e k M k\ < °-23 eV (or < 0.85 eV) at the 2a level, extracted 
from the present experimental upper bound on the Qv(3(3 decay. 

We have also looked at whether the future measurements of lepton number violation and 
non-unitarity of neutrino flavor mixing at low energies are possible to shed light on the masses 
of heavy Majorana neutrinos. Taking the minimal type-I seesaw scenario for example, we 
have illustrated the possibility of determining or constraining two heavy Majorana neutrino 
masses by using more accurate low-energy data on the Qvf3f3 decay and non-unitary neutrino 
mixing and CP violation. Such an analysis can simply be extended to the more general cases 
of the type-I seesaw mechanism with three or more heavy Majorana neutrinos. 

As stressed in Ref. [13], testing the unitarity of the light neutrino mixing matrix V in 
neutrino oscillations and searching for the signatures of heavy Majorana neutrinos N k at 
TeV-scale colliders can be complementary to each other, both qualitatively and quantita- 
tively, in order to deeply understand the intrinsic properties of Majorana particles. We 
optimistically expect that some experimental breakthrough in this aspect will pave the way 
towards the true theory of neutrino mass generation and flavor mixing. 

I am indebted to M. Chaichian for warm hospitality during my visiting stay in Helsinki, 
where this work was started, and to T. Ohlsson for warm hospitality during the Nordita 
scientific program "Astroparticle physics — A Pathfinder to New Physics" in Stockholm, 
where this work was finalized. I am also grateful to W. Chao, H. Zhang and S. Zhou for 
useful discussions and comments. This research was supported in part by the National 
Natural Science Foundation of China under grant No. 10425522 and No. 10875131. 
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FIGURES 



'Ui > tU\ 




FIG. 1. A schematic Feynman diagram for the lepton-number-violating 0vf3f3 decay, in which 
"x" stands for either light Majorana neutrinos v i (for i = 1,2, 3) or heavy Majorana neutrinos N k 
(for k = 1, ■ ■ ■ , n) in the type-I seesaw mechanism. 



